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ABSTRACT 



o 

H 

' Aims. We obtained the first long, homogenous time-series of V2I04Cyg, consisting of 679 datapoints, with the uvby/3 photometers 

' of Sierra Nevada and San Pedro Martir Observatories with the aim to detect and subsequently interpret the intrinsic frequencies of 

this previously unstudied variable star, which turned out to be a Be star We try to figure out its place among the variable B stars on 
the upper Main Sequence. In order to obtain additional information on physical parameters we collected a few spectra with the Elodie 
and FIES instruments. 

^ ' Methods. We searched for frequencies in the uvby passbands using 2 different frequency analysis methods and used the S/N>4 

CO , criterion to select the significant periodicities. We obtained an estimate of the physical parameters of the underlying B star of spectral 

. type between B5 and B7, by correcting for the presence of a circumstellar disk, using a formalism based on the strenght of the Ha 

0^ ' line emission. 

04 ' Results. We detected 3 independent frequencies with amplitudes below 0.01 mag, vi= 4.7126 d"' , V2= 2.2342 d"' and V3= 4.671 d"' , 

and discovered that V2104 Cyg is a Be star. The fast rotation (vsin i=290±10 km s"', and 27° < i < 45°) hampered the investigation of 
the associated pulsational parameters {. Nevertheless, the most plausible explanation for the observed variability of this mid-late type 
Be star is a non-radial pulsation model. 
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1. Introduction less massive B stars (spectral type B2-B9) and slow rotators that 

• pulsate in high radial order ^-mode of low spherical harmonic 

Several pulsators and variable stars of spectral type B occupy degree { with typical pulsation periods of 0.5 to 5 days. We refer 

the region of the upper Main Sequence. Among them are the /3 Stankov & Handler (2005) and De Cat (2002) for overviews 

Cephei stars, the Slowly Pulsating B stars (SPBs) and the Be j^e properties of the p Cephei stars and SPBs, respectively, 
stars. 

The former two classes are well-established pulsators that 
show multiperiodic line-profile and light variations and whose The Be stars are a more complex and less well understood 
pulsations are driven by the /^-mechanism acting on the iron class of near Main Sequence rapidly rotating B stars that pro- 
bump. Most yS Cephei stars (spectral type B0.5-B3) pulsate in duce a disk in their equatorial plane and hence show Balmer 
low-order p- and g-modes with typical pulsation periods be- line emission. Be stars exhibit variations on different time-scales 
tween 0.3 and 0.8 days. The largest fraction of this pulsational (hours-years), with a broad range of amplitudes. Explanations 
class are moderate rotators (vsin / < 150 km s '). The SPBs are variability have been sought among the dynamics of the 

circumstellar disk, stellar wind, rotational modulation, magnetic 

field and/or pulsations (e.g. Porter & Rivinius 2003). The high- 

Send offprint requests to: K. vyllsmosven -.f-.- fn.u u in. 

' est fraction of Be stars have been found among early-B type stars 
Based on observations obtained at the Observatorio Astronomico , i -i i j • . i . • . i 
. , „ r> J »j' /»j ■ N . JO- XT J (B1-B2), while the observed cases in the spectral type interval 
Nacional San Pedro Martir (Mexico), Observatorio de Sierra Nevada \. t,-, « > /r, ^^^^s f, , 
(Spain) and Observatoire de Haute Provence (France), and on observa- °^ ^^7-A2 are less frequent (Zorec 2000). However, recently a 
tions made with the Nordic Optical Telescope, operated on the island of "umber of late-type Be stars were discovered during the ground- 
La Palma jointly by Denmark, Finland, Iceland, Norway, and Sweden, based preparatory work of the CoRoT satelhte (Neiner et al. 
in the Spanish Observatorio del RoquedelosMuchachos of the Instituto 2005). The emission line B stars occupy the instability strips 
de Astrofisica de Canarias. of P Cephei and SPB pulsators and an obvious question arises 
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whether or not the observed short-period Be variations resemble 
the properties of these well-estabhshed pulsators. 

The short-periods observed in Be stars typically range from 
0.2 to 5 days. Ground-based light curves mainly reveal monope- 
riodic signals, with periods close to the rotational period (the so- 
called A Eri variables, Balona 1990), and amplitudes of a few 
tens of mmag. These targets also show line-profile variations 
(LPV) of low spherical harmonic degree, similar to the LPV 
observed in /? Cephei variables, although the latter have shorter 
periods. The variability has been attributed to rotational modu- 
lation by Balona (1990, 1995) or non-radial pulsations (NRP) 
(Baade 1982). Spectroscopically, nearly all Be stars of early-B 
type show LPV (e.g. Rivinius et al. 2003). Contrary to ground- 
based photometric data, the LPV of several early-type Be stars 
reveal multiperiodic signals. The observed variability is inter- 
preted in terms of NRP g-modes, typically with ( - m - +2. 
Given the g-mode nature and the similar length of the pulsa- 
tion periods, these Be stars seem to resemble the pulsational be- 
haviour of SPBs. Until now, no clear cases of LPV in late-Be 
type stars (later than B6) have been found (Baade 1989; Rivinius 
et al. 2003). Recently, proof of the presence of low-amplitude 
additional periods and variability of late-type Be stars has been 
given by the MOST satellite as multiperiodic signals have been 
detected in the Hght of the Be stai's HD163868 (B5Ve, Walker 
et al. 2005) and /3 CMi (B8Ve, Saio et al. 2007). This result is 
very encouraging and gives high expectations for the outcome 
of the Be programme stars of the space mission CoRoT, which 
has been launched on December 27th 2006. 

V2104Cyg (HD190397, HIP98611, a'2ooo=20''01'"45^54, 
<52000=+57°39'06'.'51, V=7.65 mag) is Hsted in the HIPPARCOS 
catalogue (Perryman et al. 1997) as an unsolved variable of spec- 
tral type AO, probably a pulsating star, located at a distance of 
approximately 617 + 263 pc. According to Grenier et al. (1999), 
the spectrum of V2104Cyg resembles a B8 Main Sequence star 
Such a late-type B star or early-type A star would be a chal- 
lenging object for many variable classes. In order to clarify this 
vague characterization we obtained a photometric time-series 
and a few spectra. From the spectra it soon became evident that 
V2104Cyg shows Balmer line emission (see Sect.lH. The study 
of V2104Cyg is particularly interesting as (multiperiodic) vari- 
ability is scarcely detected in late-type Be stars (see above). We 
present an analysis of the variable character of this newly dis- 
covered Be star and try to figure out its place among the variable 
B stars on the upper Main Sequence. 

2. Photometric observations 

V2104Cyg was observed in the framework of a double-site 
campaign in the autumn of 2005. Two twin Danish six-channel 
uvby/3 photometers were used at Sierra Nevada (SNO), Spain, 
and San Pedro Martir (SPMO), Mexico, observatories, attached 
to the 90 cm and 1.5 m telescopes, respectively. Both photome- 
ters are equipped for simultaneous measurements in uvby or 
the narrow and wide li/3 channels (Nielsen 1983). Nevertheless, 
most of the data were collected in the four uvby filters in order 
to investigate the variability behaviour of this star. Only a few 
points were obtained in H/3 at SNO for calibration purposes of 
the photometric indices. 

HR 7634 (y=6.16mag, A4Vn) and HR 7692 (y=6.19mag, 
F4V) were used as comparison stars in both sites. A total of 
414 uvby datapoints of V2104Cyg were collected on 8 nights 
at SNO, from October 3 to 26, and 265 uvby datapoints on 7 
nights at SPMO, from November 4 to 14. In total, about 54 hours 
of useful data were acquired. The light curves are given in Fig.II] 



The quality of the nights was excellent at SPMO. Bouguer's 
lines, i.e. the least-squares fit to the observations of a same 
star versus airmass for each night, yielded correlation coeffi- 
cients in the interval 0.994-1.000, whereby 23 times out of 
28 in the 0.998-1.000 interval. Mean extinction coefficients 
were 0.472 mag in u light, 0.267 mag in v light, 0.165 mag 
in b light and 0.1 16 mag in y light. Similar mean coeffi- 
cients, with standard deviation given between brackets, were ob- 
tained at SNO: 0.484 mag (0.013mag), 0.278 mag (0.009 mag), 
0.171 mag (0.007 mag) and 0.121 mag (O.OlOmag) for filters u, 
V, b and y, respectively, with the airmass ranging from about 1 .05 
to 1.5 each observing night. 

An analysis of the magnitude differences between the two 
comparison stars was carried out to get an additional feeling of 
the quality of the dataset. The overall standard deviations of the 
268 magnitude differences obtained at SPMO are: 6.7 mmag in u 
light, 2.0 mmag in v light, 1 .9 mmag in b light and 1.8 mmag in y 
light. The high scatter in the u filter is due to some leakage prob- 
lems in the photomultiplier which affected the measurements in 
a random way. The standard deviations of each night lie between 
4.6-7.5, 1.4-2.3, 1.4-2.2 and 1.2-2.1 mmag in u, v, b and y fil- 
ters, respectively. 

In case of the SNO data, we measured the two comparison 
stars 427 times and the magnitude differences yielded standard 
deviations of 6.3, 2.3, 2.3 and 2.9 mmag in u, v, b and y filters, 
respectively. Excluding the night JD 24453652 (where values are 
1.5-2.0 times greater), the standard deviations of each night lie 
between 5.5-6.7, 1.7-2.6, 1.4-2.3 and 2. 1-2.7 mmag in u, v, b 
and y filters, respectively. 

Since small differences are present in the mean values mea- 
sured at SNO an SPMO, the magnitude differences between the 
comparison stars have been aligned. The combined timeseries 
have standard deviations of 6.4, 2.2, 2.1 and 2.3 mmag in uvby, 
respectively. Therefore, our photometric timeseries are charac- 
terized by a very good precision, close to the best value we can 
obtain from ground. The corresponding periodogram is shown 
in the bottom panel of Fig. |2] 



3. Frequency analysis 

The frequency search was carried out using Lomb- 
scARGLE(Scargle 1982) as well as the least-squares power 
spectrum method (Vanfcek 1971). The latter method allows to 
detect the constituents of the light curve one by one, whereby 
only the values of the detected frequencies are introduced 
as known constituents in each new search. Such a procedure 
differs from the scargle method as it does not require any 
data prewhitening as amplitudes and phases of the known 
constituents are recalculated for each new trial frequency, 
whereby always the exact amount of signal for any detected 
term is subtracted. 

We first double-checked the reliability of the comparison 
stars. The frequency analysis of the combined timeseries be- 
tween the two comparison stars did not reveal any significant 
term. White noise is around 0.3 mmag in the vby data and ampli- 
tudes of the highest peaks are in the 0.6-0.8 mmag interval. The 
u data show a more noisy behaviour The highest peaks occur at 
different frequencies in the different filters, suggesting a random 
enhancement. As illustration, we show the frequency spectrum 
of the V data in the bottom panel of Fig. |2] Thus, we can state 
that the two comparison stars are constant at the level of a few 
0.1 mmag. 
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Fig. 1. Light curves of V2104Cyg obtained at SNO and SPMO in October-November 2005. For each night (from bottom to top) 
the u (dots), V (stars), b (crosses) and y (+) magnitudes are plotted in time (HJD-2453000). Note the difference in scale between the 
u and vby data: tickmarks in the u filter and vby filters correspond to 0.04 mag and 0.03 mag respectively. 



The frequency analysis of the V2104 Cyg timeseries turned 
out to be more complex. First, we analysed the SNO and SPMO 
datasets separately. We searched for frequencies in the interval 
0.01-25.0 d ' with a frequency step of 0.01 d We detected 
similar structures in the least-squares power spectra and the 
SCARGLE periodograms of the datasets in the different filters and 
found evidence for two frequencies centered at v\-A.l\ d"' and 
V2=2.23 d or their one-day aliases. Additional frequency 
peaks were visible but were below the significance threshold (see 
below). 

Subsequently, a preliminary solution for each dataset was 
obtained by means of these two frequencies only and the two 
datasets were re-aligned at the same mean brightness level. In 
particular, we obtained the same mean magnitude for the y data, 
within error bars, while small misalignements were measured in 
the other data. 

In a second step we analysed the aligned full dataset, consist- 
ing of 679 datapoints and spanning 42 days. Given the improved 
frequency resolution we adopted a frequency step of 0.001 d 
The spectral window of our timeseries is shown in the upper 
panel of Fig. |2] Besides vi =4.713 d"^ and V2=2.234 d"', with 
amplitudes of respectively 7.6 and 5.7 mmag in the v filter, we 
also found in all filters and using both frequency methods a third 
term V3 =4.670 d"', which has an amplitude of 2.4 mmag in the 
V time series. The amplitudes of all 3 frequencies satisfy the 
S/N>3.6 significance criterion and also the more severe S/N>4 
criterion (Breger et al. 1993; Kuschnig et al. 1997; De Cat & 
Cuypers 2003) in the vby filters. However, the low-amplitude 
of V3 does not pass S/N>4 in the noisy u filter (see last row 



Table [T). The S/N-level was computed as the average amplitude 
over a frequency interval with a width of 5 d ' in an oversampled 
SCARGLE periodogram obtained after final prewhitening. We note 
that the S/N-level is a factor 1.5 higher than the white noise de- 
tected in the comparison stars (bottom panel Fig.|2]l. Comparing 
the white noise with the residual time-series of V2 104 Cyg (last 
but one bottom panel Fig. |2]l, there might be even more fre- 
quencies present in the signal of V2 104 Cyg, with amplitudes 
currently below the significance threshold. For the moment we 
accept a model with vi, V2 and V3. We did not find a direct re- 
lation between the three frequencies, which can be considered 
independent of each other To optimise this triperiodic model 
we used a code, kindly made available by Dr Jan Cuypers 
(Royal Observatory of Belgium), that searches simultaneously 
a set of frequencies around the input values to find the multi- 
periodic combination that fits the data best. The optimised solu- 
tion is: yi= 4.7126 + 0.0004 d"', V2= 2.2342 ± 0.0005 d"' and 
V3= 4.671 ±0.001 d"^ A similar fit using a different code (mtrap, 
Carpino et al. 1987) yielded coincident results. The given ac- 
curacy of the frequencies is the average of the estimated ac- 
curacies in the uvby filters. The frequency accuracy in a par- 
ticular filter X is calculated as cry j. = '\K(i)o'sid,x/^ V(N)Ay^^AT 
(Montgomery & O'Donoghue 1999), with cTstd. v the standard de- 
viation of the final residuals, Ay^ ^- the amplitude of the frequency 
V in filter x and AT" the total timespan of the observations. The 
solution of the least-squares fit of this triperiodic model is given 
in Table[T] The phases 0, are calculated according to the formula 
A() + 2i=i,3^i(cos2:7r(f- Tojvi + (pi). Note that the phase values 
for V3 in the four passbands are in excellent agreement, support- 
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Fig. 2. Left panel: Window function (top) and scargle periodograms of V2104Cyg at different stages of the frequency search in 
the Stromgren v fiher The ampHtudes are expressed in mmag. The detected frequencies v\- A.1126 + 0.0004 d"' , V2- 2.2342 + 
0.0005 d"' and V3= 4.671 + 0.001 d"' are indicated by dashed vertical lines. The last but one bottom periodogram shows the 
final residuals after prewhitening with vi, V2 and V3. The light grey horizontal lines represent the (S/N>3.6)-level. In the bottom 
periodogram we present the variability of the comparison stars, illustrating that they are constant at the level of a few 0.1 mmag. 
Right panel: Phase diagrams of the 3 detected frequencies in the v filter Black stars correspond to observations and the grey full 
lines represent the least-squares fit. The amplitudes are expressed in mag. 



ing the reality of this small amplitude term. The triperiodic fit 
applied to the individual datasets of SNO and SPMO yielded the 
same mean brightness levels in all filters, i.e. the alignment pro- 
cedure was insensitive to V3. The scargle periodogram of the v 
filter at different stages of the frequency search and the phase di- 
agrams with the 3 frequencies are given in Fig.|2] For reasons of 
visibility of the peaks around V3 we only show the (S/N>3.6)- 
level, indicated by the light grey horizontal line. 

We also analysed the hipparcos data in search for traces of in- 
trinsic variability. The resulting scargle periodogram resembled 
very much the window function and no frequencies could be de- 
tected. The poor spectral window and the less good accuracy of 
the HIPPARCOS data did not allow to detect any of our three small 
amplitude terms. 

4. Elodie and FIES spectra: evidence for a Be star 

In addition to the photometric time-series we obtained 3 
spectra during 3 consecutive nights in June 2006 with the 
Elodie spectrograph (R~45 000) at the 193 cm telescope of the 
Observatoire de Haute Provence (OHP), France, and one spec- 



trum in November 2006 with the FIES spectrograph in medium 
resolution mode (R~45 000) at the Nordic Optical Telescope 
located at the Observatorio del Roque de los Muchachos, La 
Palma, Spain. We combined the 3 Elodie spectra, after checking 
that there are no clear differences visible, resulting in an average 
spectrum with signal-to-noise ratio of 70. The signal-to-noise ra- 
tio of the FIES spectrum is about 150. 

The Ha and H/3 line profiles clearly show a double-peaked 
emission line (see Fig.[3]l. Hence, we identify V2104Cyg as a 
Be star. The emission effect appears to be much stronger in Ha 
than in H/?. In the FIES spectra we also see a hint for emission 
in the absorption profile of Hy but this needs to be confirmed. 
A comparison between the Ha profile (upper panel Fig. [3]) ob- 
tained in June 2006 (Elodie spectra, grey) and in November 2006 
(FIES spectrum, black) shows a radial velocity shift and small 
variations in shape of the emission peaks. A dedicated spectro- 
scopical campaign is necessary in order to study the dynamical 
character of the circumstellar disk. 

As the line profiles are strongly rotationally broadened, only 
the Balmer lines. He i lines and the absorption line of Mg 11 4481 
A are prominently present and not lost in the noise of the spec- 
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Table 1. Amplitudes (in mag) and phases (in rad) of the frequencies vi= 4.7126 d"', V2= 2.2342 d"' and V3= 4.671 d"' calculated 
from the uvby time-series of V2104Cyg by means of a triperiodic least-squares fit. For each filter also the mean magnitude difference 
is given together with the r.m.s. of the residuals. The last row indicates the calculated (S/N=4)-level at V3= 4.671 d"^ Phase - 
corresponds to To = 2453647.3034. 



Term Freq. 
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gilt 
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Ampl. 

[mag] 


Phase 
[rad] 


Ampl. 
[mag] 


Phase 
[rad] 


Ampl. 
[mag] 


Phase 
[rad] 


Ampl. 

[mag] 


Phase 
[rad] 
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0.0086 


0.17 


0.0073 
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0.0071 
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±.0004 


±.0004 


±.09 
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1.4610±.0001 


1.5124±.0001 


Residual r.m.s. [mag] 


0.0065 


0.0025 


0.0026 
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S/N=4 @V3 
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Fig. 3. The Ha and HyS profiles of V2104Cyg clearly show 
emission. The black profiles are FIES profiles while the grey 
Ha profile is taken with the Elodie spectrograph. Between June 
2006 and November 2006 the emission peak shifted and changed 
slightly in shape, indications for the dynamical character of the 
circumstellar disk. 



tra. To derive an estimate of vsin / we calculated theoretical line 
profiles for a rotating star by means of the code described by 
Schrijvers & Telting (1999), whereby we assumed a neglegi- 
ble contribution to the line shape from NRP. We used limb- 
darkening coefficient a=0.444, taken from Dfaz-Cordoves et al. 
( 1 995) for the physical parameters given in Table|2] As Collins et 
al. (1991) suggest the He i 4027 A, He i 4471 A and Mgii 4481 
A profiles as most suitable candidates for the determination of 
vsin /, we focussed on these line profiles only. From the best fits 
of the observed profiles we estimated vsin/=290 km s"', with 
formal errors of 10 km s"' . 

We calculated the EW of the Ha profile and obtained EW~ 
-6.0 ± 0.5 A. The red peak of the emission is higher than the 
violet peak with peak intensities of R=2.42 and V=2.36. The 
peak separation AV between the red and violet peak is about 
105 km s"' . As AV is a measure of the radius of the circumstel- 
lar disk, namely (2 vsin //Ay)^=(/?di,sk/^*) assuming a Keplerian 



rotation (e.g. Zamarov et al. 2001), we estimate the radius of the 
circumstellar disk to be about 34 times the stellar radius. 

Given the double-peaked nature of the emission line, with 
small side-lobes on each side (Fig. [3]), we expect to see the cir- 
cumstellar disk under an inclination angle / between 45° and 
~ 20° (e.g. Hanuschik et al. 1996). 

We compared the FIES spectrum of V2104 Cyg with a set 
of spectra from the noao spectral library (Valdes et al. 2004) and 
deduced a spectral type between B5 and B7. Hence we obtain 
a hotter star than previously reported (AO, Ferryman et al. 1997; 
B8, Grenier et al. 1999). The misclassification of a late B type 
star as an early A type star is not uncommon among fast rota- 
tors. Being a mid-late Be star, V2104Cyg remains an interesting 
target given that cooler Be stars are less frequently observed than 
early-type Be stars. 

5. Physical parameters 

First, we calculated T^ff, \ogg and My of V2 104 Cyg from the 
Stromgren indices, whereby using HR 7634 and HR 7692 as cal- 
ibration stars. The observed Stromgren indices (all expressed in 
mag throughout the text) are (b - y) - 0.000, nii - 0.096 and 
C[ = 0.534. They have been calculated from the mean mag- 
nitude differences reported in Table [T] and from uvby standard 
photometry of HR 7634. A = 2.68 1 value simulaneous to uvby 
photometry has been determined from dedicated observations at 
SNO. We obtained dereddened values starting from the above 
indices and applying the TempLogG method (Kupka & Bruntt 
2001): {b - y)o = -0.057, mo = 0.115, c,, = 0.523. The inter- 
stellar reddening is non-negligible, - H-0.057mag, since 
the star lies only 14° above the galactic plane. Additionally, 
the TempLogG method allowed to derive My = -1.44 mag, 
\ogg - 3.37 and Teff = 14000 K. According to these values, 
V2 104 Cyg is located far above the instability strip of the SPBs 
in the Hertzsprung-Russell diagram (black square in Fig.lUi. 

As we are dealing with a Be star we expect that the presence 
of the circumstellar envelope affects the observed photometric 
indices. This might also explain why simbad reports scattered 
values (2.610, 2.648 and 2.705) of the /5 index for our target. To 
deduce the indices from the underlying B star we calculated the 
transformation formulae given by Fabregat & Reglero (1990), 
which are based on the assumption that the contribution of the 
circumstellar disk is proportional to the EW of the Ha emission 
peak. Assuming an EW of Ha of -6.0 + 0.5 A we obtained 



6 



K. Uytterhoeven et al.: Multiperiodicity in the newly discovered mid-late Be star V2104Cygni 



Fig. 4. Possible position of V2104Cyg in the Hertzsprung- 
Russell diagram. The black square indicates the position derived 
from the observed photometric indices while the star indicates 
the position after correction for the contribution of the stellar 
disk. The full black lines and the dashed black lines represent the 
instability strips for the SPB and yS Cephei modes respectively, 
with a frequency between 0.25 and 25 d"' and { < 4 computed 
for Main Sequence models (De Cat et al. 2007). The dark grey 
dotted lines show the ZAMS and the TAMS. 




the following new values of the Stromgren indices: (b - y)* - 
-0.01 ±0.01, (mi)* = 0.10 ±0.05, (ci)* = 0.57 + 0.05 and = 
2.73 + 0.01. Applying the TempLogG method on these values 
resulted in dereddened values (/? -y)(),* - -0.05 ± 0.01, (mo)* - 
0.12+0.05 and (co)* = 0.564 in addition to (M^)* = -0.49 mag, 
(logg)* = 3.94 ± 0.05 and (^eff)* ~ 13 500 K. After correction 
for the contribution of the circumstellar disk, the position of the 
underlying B star falls nicely inside the SPB instability strip (★ 
in Fig.m. 

From the corrected value My,* = -0.49 mag we can estimate 
the luminosity and the mass of V2104Cyg using the relation for 
Main Sequence stars between absolute magnitude and luminos- 
ity and the mass-luminosity relation. Assuming that V2104Cyg 
is a MS star we find: L/Lq ~ 127 and M/Mq ~ 4. Using the 
formula log{R/Ro) = ^\og(L/LQ) - llogiT^fi/T^fi^Q), we esti- 
mate R/Rq ~ 2.1. The value of the mass is compatible with the 
expected values for a B5-B7 star according to the tabulation of 
Harmanec (1988), while we obtained a slightly lower radius (2.1 
versus 3.0 /?©). From these estimates of mass and radius, the ex- 
pected upper limit of the inclination ; < 45°(see above), and with 
vsin /=290±10 km s"', we can derive that V2104Cyg rotates at 
least at 65% of its critical velocity, using Vcrit - y/GM^jR^. This 
fact confirms the tight connection between fast rotation and the 
Be phenomenon. Moreover, an interpretation in terms of mul- 
tiperiodic pulsations would confirm and strenghten the relation 
between pulsation and mass loss, and demonstrate that this re- 
lation also exists in mid-late Be stars. Subsequently, we obtain 
from the derived vsin / interval and imposing v<Vcrit a lower limit 
for the inclination angle: / > 27°. 

Next, from vsin /=290±10 km s"', a stellar radius in the in- 
terval [2,3] Rq, and an inclination angle between 27°and 45° 
we arrive at a (very) rough estimate of the rotational frequency: 
Q 6 [2.6, 6.5] d"'. We notice that the observed intrinsic frequen- 
cies vi and vj lie in this interval. It is not uncommon for Be 
stars that the rotational period is dominantly present in the light 
variations (e.g. Balona 1990). We will discuss this possibility 
further in Sect. |6] 



Table 2. Estimates of the physical parameters of the Be star 
V2104Cyg, after correction for the contribution of the circum- 
stellar disk. The eff'ective temperature is given in K; the mass 
and radius are expressed in solar units; vsin / is given in km s ' . 



' cff 



Mass 



Radius 



13 500 ± 1350 3.94 ±0.39 ~4 2.5 ±0.5 



127 290±10 



We stress that the derived physical parameters tabulated in 
Table|2]are only a first estimate as the contribution and properties 
of the circumstellar disk are still unknown and as the properties 
derived from the low S/N noise spectra have to be taken with 
caution. 



6. Discussion and conclusions 

In the light curves of V2104Cyg, we detected 3 intrinsic fre- 
quencies with amplitudes below 0.01 mag: vi= 4.7126 d"', 
V2= 2.2342 d"' and V3= 4.671 d"'. We also discovered that this 
object is an intrinsically fast rotating Be star V2104Cyg turned 
out to be a particularly interesting case as there are not many 
precedents for detection of multiperiodicity in mid-late type Be 
stars based on ground-based photometric observations. Indeed, 
from previous efforts (Cuypers et al. 1989; Balona et al. 1992) it 
appears to be not obvious to detect multi-periodic signals from 
ground-based photometry in Be stars. Only recently Gutierrez- 
Soto et al. (2006) claimed to have detected unambiguously sev- 
eral NRP periods in the light curves of NW Ser (spectral type 
B3) and V1446 Aql (spectral type B5). If V2104Cyg is a multi- 
periodic non-radial pulsator as well it would be the third known 
Be pulsator with multiperiodic signals detected in its ground- 
based light curves. 

How can we explain the variability of V2104Cyg? Short pe- 
riods of 0.2"', like v\ and are typically detected in early-B 
type stars such as (3 Cephei pulsators or early-Be stars, but are 
quite unusual in late-type B(e) stars. As V2104 Cyg lies in the in- 
stability domain of the SPBs, far from the overlap domain with 
the j6 Cephei pulsators, we expect to detect longer periods, such 
as V2- Rivinius et al. (2003) identified {-m-2 g-modes associ- 
ated to periods longer than 0.4"' in a selected set of Be stars. 
Given that V2 104 Cyg has a spectral type between B5 and B7 it 
is worth noting that the NRP stars in Rivinius et al. (2003) have 
a spectral type earlier than B6. 

According to Saio et al. (2007) it is likely that all rapidly 
rotating Be stars show NRP. For fast rotating stars, the com- 
parison of the observed frequencies Vobs with the eigenfrequen- 
cies of the star Veigen provided by theoretical calculations is 
not straightforward (e.g. Berthomieu 1978; Dintrans & Rieutord 
2000). The rotation lifts the degeneracy of the eigenfrequencies 
of the modes in m and gives rise to a perturbation of the fre- 
quencies (e.g. Saio 1981; Dziembowski & Goode 1992; Soufi 
et al. 1998). These frequency shifts depend on w, the rotation 
frequency Q, and other parameters, and can be very large, es- 
pecially for g-modes. This perturbation can explain why the ob- 
served frequencies v\ and V3 of V2104Cyg are larger than ex- 
pected for stars within the instability strip of SPBs. Moreover, 
when the oscillation frequency is of the same order as the ro- 
tational frequency, the whole perturbative approach is invalid 
and the familiar structure of modes (n, {) disappears (Dintrans 
& Rieutord 2000). Hence, rotation has a significant effect on the 
non-adiabatic analysis and mode identification. Consequently it 
is not easy to assign /-values to the observed frequencies of a 
fast-rotating star, like V2104Cyg. In particular, the fast rotation 
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of V2104Cyg prevents the application of a mode-identification 
method that does not account for higher order rotational effects, 
such as the method of photometric amplitude ratios made avail- 
able by Dupret et al. (2003). A discussion of the implications of 
applying standard mode-identification techniques, used to iden- 
tify modes in non-rotating stars, to fast-rotating stars can be 
found in Townsend (2003). 

It is interesting to point out that for V2104Cyg the ampli- 
tudes associated to V3 have, unlike the ones associated to vi and 
V2, a similar value in all four filters (see Table[T]). In general, the 
amplitudes of NRP modes in B stars decrease from the u filter 
towards the vby filters. However, similar pulsation amplitudes 
at different wavelenghts are predicted for some of the modes of 
theoretical models. Several examples can be found in De Cat et 
al. (2007; e.g. HD 14053 and HD 89688). 

An alternative explanation of the observed variability of 
V2104Cyg might be in terms of a model based on rotational 
modulation. In many Be stars the detected period is consistent 
with the expected period of rotation of the star (Cuypers et al. 
1989; Balona et al. 1992). As seen in Sect.|5]might vi or V3 be 
related with Q. Our dataset did not allow a tight enough con- 
straint on O to confirm this. The frequency V2, however, is much 
lower than the expected rotational frequency. Therefore, a model 
based on rotational modulation could not explain the light vari- 
ability of V2104Cyg. 

The unambiguous detection of photospheric multiperiodic- 
ity in a Be star is in favour of a NRP interpretation (e.g. Porter 
& Rivinius 2003). As we detected 3 independent frequencies in 
the light curves of V2104Cyg, the NRP hypothesis seems cur- 
rently the most plausible explanation. The interpretation of pul- 
sational variability in Be stars is complicated by fast rotation. 
Recently, theoretical progress on this topic has been made by 
Townsend (2005), who proposed a new explanation for unstable 
modes in mid- to late B-type SPB and Be stars in terms of retro- 
grade mixed modes instead of ^-modes and by Dziembowski et 
al. (2007), who found theoretical backing for the observed fre- 
quencies in the MOST data of the Be stai- HD163868 (Walker et 
al. 2005). The case study of V2104Cyg, a fast-rotating mid-late 
Be star with 3 excited modes, is another challenge for theoreti- 
cians in the game of matching theory with observations. 

A high-resolution spectroscopic time-series of V2104Cyg 
could be a future, possible observational progress. As LPV 
are rarely detected in late-type Be stars (Baade 1989; Rivinius 
et al. 2003) it would be an interesting excercise to see if 
V2104Cyg is an exceptional case with this respect and indeed 
shows LPV. If present, a line-profile analysis and subsequent 
mode-identification will allow to assign m-values to the detected 
modes. Additionally, a study of the behaviour of the Ha emis- 
sion peak will allow a study of properties and behaviour of the 
circumstellar disk. 
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